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In this talk I present the current status of the family of Monte Carlo programs for 
four-fermion physics based on the Yennie-Frautschie-Suura resummation of soft 
real and virtual photons. I focus on their applications to LEP2. 



1 Introduction 

The role of Monte Carlo (MC) codes at LEP2 is more important than at 
LEP1. This is primarily because the physics is more complicated: we deal 
with a multitude of four-fermion final states, their topology is complicated and 
not that many semi-analytical results exist. As a consequence, for example, 
the W-mass fits are performed at LEP2 with the help of MC. At the time of 
the LEP2 Workshop in 1995 a list of wishes concerning the features of the 
"ultimate Monte Carlo" was constructed, and the goal of 0.5% precision tag 
for the total cross-section for WW-physics was set El. As the experiment is 
already well advanced it is time to confront these expectations with reality. I 
will do it for the family of four-fermion MC codes based on the YFS technique: 
KoralW, YFSWW and YFSZZ developed by our group. Let me begin with a 
small technical introduction. 
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2 YFS Monte Carlo Approach to Four-Fermion Calculations 

In their classical paper i, Yennie, Frautschie and Suura (YFS) reorganized 
the perturbative series of QED photonic corections for an arbitrary process 
in a manifestly infrared (IR) finite way. The key step was in the extraction 
of the universal real (S) and virtual (S) functions containing IR singularities. 
After properly integrating and resumming them to all orders the remaining new 
perturbative series in n functions became IR-finite. As compared to the fixed- 
order approach this layout is especially convenient for MC algorithms. One 
is not faced here with problems of IR cut-off, negative distributions or large 
real-virtual cancellations. Moreover, one can generate a multiple photonic 
bremsstrahlung instead of just one or two photons at hiqsjL, This approach has 
been pioneered by two of us (SJ and BFLW) in refs. EjulI. In mathematical 
terms, thfij-master formula of YFS series, adapted to the MC needs, is the 
followingLrl 
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To explain it very briefly, the exp(. . . ) function is the place where IR virtual 
C3tB) and real (/ S) singularities cancel; Y[d 3 kiS describe the multiple pho- 
tonic emission and ^ Pi is the perturbative expansion of a non-IR part of the 
matrix element (truncated to the third order). To lowest order, (3q is just 
the Born matrix element. This formula is the basis of the MC codes that we 
developed for four-fermion physics. 

3 The "Four-Fermion MC Toolbox" 

For the moment we provide three complementary MC codes: KoralW, YF- 
SWW and YFSZZ. I will refer to them together as "Four-Fermion MC Tool- 
box" . 

• KoralW 

The main features of KoralW BLrQ are: (1) it generates all four-fermion 
final states with the complete massive Born-level matrix elements and two 
pre-samplers for complete, massive, four-fermion phase space. The matrix el- 
ement comes from the automated package GRACE v. 2.0. (2) Anomalous 
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WWV couplings are included in CC03 graphs. (3) Multiple initial-state pho- 
tons with finite px are generated by the YFS technique and a QED 0(a 3 ) 
leading-log matrix element is included. (4) Coulomb effect," Naive QCD" cor- 
rection and non-diagonal CKM matrix are included. (5) JETSET, PHOTOS 
and TAUOLA are interfaced. (6) Semi-analytical routine KorWan for CC03 
graphs is included. (7lAnalysis of Bose-Einstein effect is done as a stand-alone 
application (in C++) El. 

As compared with the previous version (1.33) the main novelties of v. 
1.41 are the following: (1) we have added a second, independent presamplcr 
for the four-fermion phase space, which becomes a very powerful test of the 
code; (2) we have added the third-order leading-log correction to the QED 
ISR matrix element; (3) the anomalous couplings can now be parametrized in 
three different ways; (4) any combination of final states can now be specified by 
the user; (5) CKM matrix and colour reconnection probability are now input 
parameters. 

The most difficult part of the code is the phase-space generator. The 
origin of the difficulty is in the multitude and complexity of the available final 
states. Naively counting one has to deal with 9x9 decay channels of WW 
(CC processes) and 11 x 11 of ZZ (NC processes), 202 in total. Each channel 
can contribute up to 100 Feynman graphs, so the total number of "objects" to 
integrate exceeds 10 000. The solution lies in a multi-branch MC algorithm. 
This approach can already be found in TAUOLABeI and FEPlMISVEI MC 
codes. The general formula is the following: 
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Skipping details (seei), let me only say that the process specific information is 
hidden in the p CR functions. These functions must contain all mass and angu- 
lar singularities of the Feynman diagrams: 1/s, l/[(s — M 2 ) 2 + M 2 T 2 ], 1/i, 1/u 
and so on. All the above simple functions are used in f l CR in the form of sub- 
sequent sub-branches with their own splitting probabilities. It is a non-trivial 
task to fine-tune all these coefficients to get the efficient modelling of all the 
matrix elements. 
• YFS WW 

Starting from the same master formula, the YFSWW codel3oE3 features 
the O(a) corrections to W-pair production: (1) CC03 (signal) graphs are in- 
cluded in the matrix element; (2) first-order corrrections for WW production 
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process are included in the_exact form of refs.E^O as well as in an improved 
Born approximation of ref. IB; (3) YFS-type bremsstrahlung is generated from 
both initial and intermediate (WW) states; (4) anomalous couplings are in- 
cluded in the Born matrix element; (5) JETSET, PHOTOS and TAUOLA are 
interfaced. 

On the technical side, YFSWW presents a number of new ideas. First 
of all the YFS scheme originally derived for fermions had to be extended to 
bosons. In a nutshell this is possible since soft photons are "blind" to spin. 
Both the fermionic and bosonic vertices look the same in the IR limit: 
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The second important extension of the original YFS paper □ was in deriving 
the YFS form factors in the case of heavii. massive particles, contrary to the 

Siginal small mass approximation of ref. □. These functions can be found in 
. Next, since photons are radiated from the W-bosons with a finite width, 
one must do it in a way that respects gauge invariance. We did it by adding 
compensating, loop corrections that restore gauge invariance in a similar way 
as in refsliSEi Finally, we had to avoid double counting of Coulomb effect in 
the matrix element and in YFS virtual B-function, both arising from the same 
type of loop corrections. We have done that by a proper redefinition of the 
B-function. 

With the help of the YFSWW code we can evaluate the size of 0{a) 
corrections, see the table below (on the example of the csev e final state). 



Ecu [GeV] 


cr [pb] 


(<rf - a^)/a 




161 


0.1768 


-0.83% 


+0.22% 


175 


0.5891 


-1.32% 


-0.006% 


190 


0.6792 


-1.71% 


-0.22% 


205 


0.6850 


-2.22% 


-0.61% 


500 


0.2710 


-4.72% 


-3.09% 



The main novelty of this result is that the corrections are calculated within 
the realistic full four-fermion MC framework, allowing for any experimental 
cuts. The numbers show that the 0(a) correction, of about 2%, is similar 
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in size to the e + e~ — > W + W~ case analysed in the literature, seetl and refs. 
therein. Therefore it must be included in the MC in order to reach the targeted 
0.5% precision level. The approximate Born cross-section barely fits within the 
0.5% limit, depending on the actual highest energy available at LEP2. 

The last question to be addressed here is how to combine various phys- 
ical corrections of KopJW and YFSWW? This almost trivial question re- 
quires some attention c3 since certain effects, if they are put together im- 
properly, can lead to .severe effects, such as the case of running W-width 
in background graphs E3. The natural way out is to use a "Multiparameter 
Linear Interpolation" , that is just a series expansion around the "minimal" 
CC03 cross-section: a M Li = CCC03 + (cr C caii ~ OCC03) + (o"o(a) ~ vccoz) + 
. . . (r(s), ACC, NQCD, . . . ). As the above expression has quite a few terms it 
becomes impractical in actual usage. We have shown in ref . c3 an example of 
the total cross section (with some cuts) that one can equivalently use a "Com- 
mon Sense Interpolation" by switching on all available corrections in the same 
MC run, provided one does some minimal cross-checks for each type of ob- 
servables. This way the above expresion would be reduced to two corrections 
only (from KoralW and YFSWW), and one day, if also 0(a) corrections are 
included in KoralW, to just one MC run. 

• YFSZZ 

YFSZZ E3 is a dedicated code for the production and decay of Z-pairs in 
the process e+e" -> ZZ -> hfxhh- lt features: (I) signal NC02 (ZZ doubly 
resonant) matrix element; (2) anomalous ZZV couplings; (3) multiple YFS- 
based bremsstrahlung with an 0(a 2 ) leading- log QED matrix element. This 
code is yet to be explored and developed as the physics of Z-pairs progresses. 

4 Precision of "Four-Fermion Toolbox" 

We can now turn to the basic question: What is the overall precision of the 
(KoralW+ YFSWW) "Toolbox" for the WW physics at LEP2? As stated ear- 
lier, we focus on the total cross-section with an eye on a 0.5% precision level. 

• KoralW 

(1) Technical precision is estimated at 0.2% based on: (a) internal compar- 
isons of two presamplers, (b) comparisons with other codesEJEiJ, (c) comparisons 
with semi-analytical results of the KorWan code (for CC03 matrix element). 

(2) Physical precision is estimated at 2% based on the size of the O(a) 
correction calculated by the YFSWW code. 

• YFSWW 

(I) Technical precision is estimated at 0.2% by comparing results of two 
technically different implementations of the code: YFSWW-2 and YFSWW-3 
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as well as comparison with KoralW. 

(2) Physical precision is still under study. Preliminary estimate of 0.5% 
uses in particular the fact that the so-called "non-factorizable" corpaticps are 
negligible. This was first pointed out by Fadia,,, Khoze and MartincaJHSEJ and 
then worked out in detail by other groups c3'E3'E§E3l2j. 

• Total 

Collecting the above uncertainties, we can give the preliminary estimate 
of the precision of the (KoralW + YFSWW) "Toolbox" for the total cross- 
sections of the WW-physics at LEP2 to be 0.5% (preliminary). 



5 Under the carpet 

From previous sections one may get the impression that our "Four-Fcrmion 
MC Toolbox" is almost finished. Unfortunately this is not true, especially for 
the NC processes. I will present some of the outstanding problems. They are of 
a universal type and can be even several times bigger than the 0.5% precision 
tag, see ref. a for more comments. 
• Numerical instabilities 

The ratio m|/ 's is at LEP2 of the order of 10 -12 . Together with delicate 
gauge and unitarity cancellations it leads to numerical instability problems in 
the matrix element calculations. Let me give an explicit example with an event 
generated by KoralW in the e~D e v e e + channel: 

pdg p_x p_y p_z E 

11 -.000000341278492 -.000001614768132 92.983165682216736 92.983165683620868 
-12 -.633000329710363 .109862863634447 -11.985324403772751 12.002531455067720 

12 .771573326908979 .358415826559387 -70.475962939695890 70.481097753801436 
-11 -.138572655920124 -.468277075425702 -10.521878338748095 10.533205107509984 

four-fermion weight = 913570469940928.500 

The outgoing electron is highly collinear, the corresponding transfer small, 
and the four-fermion weight (i.e. matrix element) is huge. Now, let us modify 
by hand the last two digits of the p z components of four-momenta and rerun 
the event: 

pdg p_x p_y p_z E 

11 -.000000341278492 -.000001614768132 92.983165682216722 92.983165683620868 
-12 -.633000329710363 .109862863634447 -11.985324403772731 12.002531455067720 

12 .771573326908979 .358415826559387 -70.475962939695876 70.481097753801436 
-11 -.138572655920124 -.468277075425702 -10.521878338748115 10.533205107509984 

four-fermion weight = 25094.3831593953582 

The four-fermion weight (matrix element) has changed by 11 orders of 
magnitude! What can we do about this? We can use a quadruple precision. 
Our experience shows that it cures the problems. However it requires a complex 
quadruple precision that is not available on all platforms. Also the speed of the 
calculation is much lower. For that and other reasons, specified later on, we 
prefer another solution - impose additional post- generation cut-offs on these 
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dangerous corners of the phase space: 

(1) In the case of the CC-type processes it is only a tiny angular cone 
around the beams, of the size of 10 -6 rad. It influences the total cross-section 
below 0.2% and can be hidden in the technical precision for the time being. 

(2) In the case of NC-type processes the situation is much worse. There are 
two regions (for final states with the e + e _ pair) that we had to cut out: when 
invariant mass of produced pairs is smaller than v8 GeV or when the sum 
of transverse momenta squared of visible particles is smaller than 600 GeV 2 . 
Note that this latter cut allows an e + e~ pair to go to the beam pipe. 

• Limitations of ISR 

The ISR-based description of bremsstrahlung breaks down for final states 
with at least one electron (positron) collinear to the beam. For such processes 
the t-channcl photon exchange is dominant. In other words the bremsstrahlung 
is governed by log^l/m 2 instead of log s/m 2 . As a result of the missing 
log \t\/ s, too much radiation is generated, especially in the high-p^ region. De- 
pending on the actual cuts this can significantly affect the cross-section. How 
can it be cured? First of all by adding all the radiation (interferences) missing 
in the soft limit. This means that the YFS program should be applied to all six 
external particles. Although radiation in the two s- and i-channels separately 
is already implemented in MC codes (YFS3 and BHLUMI, respectively), it is 
nonetheless not easy technically to merge them into one coherent algorithm 
in the case of four-fermion final states. (In. the case of two fermions, this has 
already been done in BHWIDE MC codeEd by some of us.) 

However, even if we have the radiation corrected in the soft limit we would 
still be missing part of the hard photonic corrections. Can these be dangerous? 
Consider the following event in the LAB frame (e + e~dd from KoralW) with 
a hard transverse photon (note that this event passes the cuts defined in the 
previous section): 

================ LAB frame =================== 

pdg p_x p_y p_z E 

PHO 5.91545655722838 47.82076396799343 2.67297019944074 48.25932927890828 
PH0 -.00000027637796 -.00000037152652 -.14955853849146 .14955853849218 
PHO .01336462801524 .00882331258921 .02409193017994 .02892896864438 
PHO .00000005382246 .00000005171260 .00563331053129 .00563331053178 
1 -.24118785025702 -1.90677419291787 -2.58051680257066 3.21762744428852 
-1 -.01773699286484 -.01532459077309 .98308780906072 .98341806182419 
11 -.11234211241339 .42661722421839 60.88292736560235 60.88452568487788 
-11 -5.55755400715285 -46.33410540129616 -61.83863527375292 77.47097871243272 
angles of decay products with resp to beams: 
d quark d~ quark e- e+ 

-.80199753074017 .99971586568334 .99997374838494 -.79821678134001 

The problem appears when we transform it to the "effective frame" . In 
order to calculate the matrix element that is defined just for four- fermions, 
the momentum carried out by photons has to be compensated for and some 
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"effective" on-shell eiectron "beams" have to be constructed. After the trans- 
formation the event becomes: 

============ effective CMS frame ============= 

pdg p_x p_y p_z E 

1 -.008731476493231 .021808545200827 -2.679767691308820 2.679889313144880 
-1 -.016063888417215 -.021968761695611 1.057202898562928 1.057600417800303 
11 .000815509213020 -.001379443999798 66.011587576947306 66.011587598375669 
-11 .023979855697426 .001539660494582 -64.389022784201416 64.389027269943767 
angles of decay products with resp to beams: 
d quark d~ quark e- e+ 

-.99996157863668 .99966881915163 .99999999970535 -.99999993036524 

As we see, we have created a "monster" - a very "small angle" configu- 
ration out of a quite transversal event. This leads to a physically unjustified 
enhancement in the matrix element calculation, i.e. to high weights or even nu- 
merically unstable ones. What is the solution? Surely, one needs the complete 
0(a) or even higher order photonic corrections. This is a very difficult task 
for the future. For now we use temporary fix-up by cutting out the high px 
photons with X) Pt — 300 GeV 2 . Unfortunately, it is not a completely effective 
cut and on an occasion numerical instabilities survive. 



6 Summary 

In this talk I presented the "Four-Fermion MC Toolbox" . It consists of three 
MC programs: (1) KoralW with all four-fermion processes, YFS-based ISR 
and anomalous WWV couplings; (2) YFSWW with CC03 matrix element, 
O(a) electroweak corrections to W-pair production, (Initial+WW)-state YFS- 
based brcmsstrahlung and anomalous WWV couplings; (3) YFSZZ with NC02 
matrix element, anomalous ZZV couplings and YFS-based ISR. The current 
precision of the codes for the WW physics at LEP2 is 2% for KoralW and 
0.5% (preliminary) for YFSWW for the total cross section. Based on this 
we expect the corresponding total precision of the Toolbox to reach the gold- 
plated 0.5% level soon. The important limitations that need to be resolved in 
the future concern the final states with electrons (positrons) at small angles, 
that is mostly the NC-type processes. They include: (1) numerical instabilities 
in matrix element (under control by applying cuts or by the use of a quadruple 
precision); (2) lack of "i-channel" bremsstrahlung; (3) lack of hard-photon 
matrix element. We are on the way to resolving these problems, and we are 
looking forward to it with excitement. 
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